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Abstract 

The electrical resistivity of a pure sample of a thin metallic film is found to 
depend on the boundary conditions. This conclusion is supported by a free- 
electron model calculation and confirmed by an ah initio relativistic Korringa- 
Kohn-Rostoker computation. The low-temperature resistivity is found to be 
zero for a free-standing film (reflecting boundary conditions) but nonzero 
when the film is sandwiched between two semi-infinite samples of the same 
material (outgoing boundary conditions). In the latter case, this resistivity 
scales inversely with the number of monolayers and is due to the background 
diffusive scattering by a finite lattice. 
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I. INTRODUCTION 



It is well known that the low-temperature electrical resistivity of an infinite pure metal 
is essentially zero because the Bloch waves associated with the underlying periodic structure 
undergo no diffusive scattering in the absence of impurities and imperfections. Then, one 
could ask what the corresponding result would be for a thin metallic film. This question, 
which has not been addressed in the literature, is the central problem tackled in this paper. 
Our conclusion is that the electrical resistivity of a thin film depends upon how the exper- 
imental situation is set up, namely, upon the boundary conditions; in particular, with an 
appropriate choice of boundary conditions, the resistivity is not zero. 

The problem of the electrical resistivity of a pure sample of a thin metallic film is of 
great current interest. It can be viewed as a fundamental question for which an exploratory 
analytic calculation may provide deeper insight as well as a limiting test of complex ab 
initio methods. In fact, this question can be posed in the context of the first-principle fully 
relativistic layered version of the Korringa-Kohn-Rostoker (KKR) method, which has been 
recently developed 0, to perform ab initio computations of magnetotransport in magnetic 
multilayers in the coherent potential approximation. In effect, this method can be easily 
applied to the computation of the electrical resistivity of a thin film consisting of a finite 
number M of monolayers of a metal, by properly specifying the conditions enforced at the 
boundaries of the system. The first result was obtained for a thin film of copper sandwiched 
between two semi-infinite metallic blocks of the same material; the resistivity was computed 
for a film consisting of between M = 6 and 45 monolayers and for current in the plane of the 
layers (CIP) and was found to decrease approximately as 116/M /zficm, for M > 20 ML — a 
result that clearly extrapolates to zero only as M approaches infinity. This nonzero value is 
somewhat unusual and led us to ask that the computation be done under different boundary 
conditions. For a free-standing slab, namely, when the copper film is surrounded by vacuum 
on both sides, it was found that the resistivity is essentially zero 0. In conclusion, when 
one performs these two computations with unequal boundary conditions, one gets unequal 
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resistivities. In other words, the resistivity depends on the boundary conditions. 

The results described above raise a number of questions that we address in this paper. 
Our goal is to explore electrical conduction in a pure sample of a thin metallic film and resolve 
the following fundamental issues, (i) Is the resistivity indeed dependent on the boundary 
conditions? (ii) How can one understand the finite resistivity of a thin, yet otherwise perfect, 
film? (iii) How can one measure the finite resistivity of a thin film? 

As we will see in Sec. |T[ a proper understanding of boundary conditions is a prerequisite 
for a thorough discussion of these questions. This analysis will be followed in Sec. [Hi] by the 
application of the Kubo formula to the case with outgoing boundary conditions, in Sec. [TV] 
by the mathematical characterization of the finite periodicity of the lattice, in Sec. [V] by 
the implications of the film finiteness on transport properties, and in Sec. |V| by concluding 
remarks. 

II. RESISTIVITY AND BOUNDARY CONDITIONS 

Electrical conduction in a thin metallic film can be modeled by representing the film in 
terms of a finite number M of monolayers or perfect atomic planes arranged periodically 
in the direction perpendicular to the boundaries of the film. Moreover, a pure sample is 
characterized by the absence of impurities or imperfections; in addition, at low temperatures, 
other resistivity sources are rendered ineffective. 

At first sight, one might just apply the standard folklore, namely, that quantum- 
mechanical Bloch waves in a periodic structure propagate without electrical resistance. In 
effect, electrical resistance arises from the loss of linear momentum information due to dif- 
fusive scattering, i.e., scattering that randomizes the electron's momentum, so that the 
outgoing electron has no knowledge of the direction of its incoming momentum [Q. In fact, 
in a pure metallic sample, electrons undergo Bragg scattering, which being highly direc- 
tional, is ineffective as a momentum-randomizing mechanism (as we will discuss in greater 
detail in Sec. |Vp. This would seem to imply that the resistivity of a perfect film should be 
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identically zero; however, this line of reasoning is simplistic: the standard folklore applies 
only to an infinite periodic sample, for which the electrical resistivity is indeed zero, un- 
der the conditions described above. Instead, when one analyzes a finite film, the question 
"What is the resistivity of a thin metallic film?" is immediately replaced by "Is the film 
really periodic?" Then, from a mathematical viewpoint, the film fails to be strictly periodic 
because of the boundaries; by abuse of language, one could describe the film as exhibiting 
"finite periodicity." It turns out that the finiteness of the film implies the existence of back- 
ground diffusive scattering in addition to ordinary Bragg scattering (see Sec. |V|), and it is 
this diffusive scattering that becomes the source of a finite-size resistivity if the boundary 
conditions are appropriately selected. 

The simplest boundary condition is provided by the ideal free-standing slab when a 
film with perfect boundaries is inserted in a vacuum. Then, the electrons undergo specular 
reflections at the boundaries (with infinitely high potentials representing the onset of a 
vacuum) and effectively "probe" a truly periodic potential. This amounts to repeating 
the film periodically; periodicity is restored by the boundary conditions and the standard 
folklore applies: the resistivity p is indeed zero. In fact, the absence of electrical resistivity 
can be traced back again to the electrons keeping their memory of linear momentum. In 
the old Fuchs-Sondheimer transport model [[|, this reflecting boundary condition, which 
amounts to the absence of diffuse scattering at the boundaries (100% specular reflection), is 
parametrized by p = 1, where p is the coefficient characterizing the specularity of scattering 
off the surfaces of the film. The ensuing resistivity p = has been confirmed by the 
relativistic layered KKR method 0. 

Does this mean that the resistivity is always zero? Of course not; the same relativistic 
layered KKR method showed that finite-size effects are not negligible when the film is 
sandwiched between two semi-infinite samples of the same material. The novelty lies here in 
the use of different boundary conditions. In effect, if the boundary conditions dictate that 
the electrons cannot keep their memory of momentum, then the background scattering of 
the finite lattice becomes the source of a nonzero electrical resistivity. 
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An extreme form of this loss of momentum information is achieved when the momentum 
of electrons entering the finite sample is totally uncorrelated with that of electrons leaving 
the sample. One can conceptualize this extreme loss of momentum information in three 
equivalent ways: in terms of reservoirs, in terms of boundary scattering, and in terms of 
boundary conditions for the electron propagators (Green's functions). First, momentum loss 
is effectively implemented by having reservoirs that absorb the electrons upon leaving the 
sample — when the electrons probe a reservoir, they get out of synch with respect to their 
"proper" behavior in the sample. Second, from the scattering viewpoint, this information 
loss can be modeled by perfectly diffuse scattering at the boundaries of the film; this con- 
dition is precisely equivalent to the choice p = in the Fuchs-Sondheimer model 0] (0% 
specular reflection, a condition that totally erases momentum memory). The third view- 
point is needed when applying the Kubo formula in terms of electron propagators; in our 
original theory of transport in metallic superlattices ||, we described the required condition 
associated with this momentum loss as outgoing boundary conditions. 

Transport theory with outgoing boundary conditions is based on the following ideas. 
Transport is described via retarded Green's functions that represent the propagation of 
electrons in the environment provided by the sample, with the condition that the electrons 
leave the film irreversibly at the boundaries. This implies the use of Green's functions cor- 
responding to an infinite medium (bulk) — at this level the calculation does not acknowledge 
the finiteness of the film whose conductivity is calculated. Then, as the propagators them- 
selves do not satisfy boundary conditions that keep track of the momentum of electrons 
scattered within the film, the calculation yields a nonzero resistivity (see Sec. [V|). In other 
words, all finite conductors have self-energy terms in their propagators that describe their 
contact with leads or reservoirs. The resistance, which is proportional to the imaginary part 
of the self-energy, reflects the fact that an electron in a finite conductor will eventually leak 
out into the leads attached to it ||. 

The discussion above has dealt successfully with the first two questions posed in Sec. |: 
the resistivity is indeed dependent on the boundary conditions and we have understood 
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conceptually how the finite resistivity of a perfect thin film arises under outgoing boundary 
conditions. However, the issue of how to measure the finite resistivity of a thin film has 
not yet been clarified. From the experimental viewpoint, the finite electrical resistivity of 
a perfect thin film still remains puzzling, even when the concept of finite periodicity is 
introduced. In effect, if we isolate the film from a bulk sample of copper and maintain it in 
contact with the "remainder" of the system (two perfectly conducting semi-infinite copper 
blocks), then, the current is shunted by these contacts. In other words, it would seem that 
it is not possible to measure the resistivity of the film in this straightforward way. In fact, 
one concludes that, for a measurement of electrical resistivity, the boundary contacts should 
not be of the same material. 

So how does one observe the calculated resistance? As mentioned in the previous para- 
graph, the shunting of the current by the contacts prohibits one from measuring the resis- 
tance of a portion of a metal. However, if the current probe is sufficiently narrow to contact 
only the finite layer, and appropriate boundary conditions are enforced, then the current 
will be limited only to the finite layer. Then, one can either measure the resistance of a 
free-standing film or have it supported on an insulating substrate. If the boundaries are 
ideal so that they have no appreciable roughness, i.e., for p — 1, the boundaries simulate 
reflecting or free-standing boundary conditions and lead to zero resistivity. Alternatively, if 
the boundaries are sufficiently roughened, they simulate current flow only in the finite layer 
subject to the p = boundary condition — this amounts to outgoing boundary conditions; in 
fact, the roughened interface could even separate the finite layer from a semi-infinite sample 
of the same material on either side, a situation that is modeled by the corresponding rela- 
tivistic KKR computation ||. It is in this way that we understand the paradox of using a 
film with a rough surface to measure the resistance calculated for a perfectly flat film; it is 
the boundary condition for transport on the surface of the film that happens to be the same 
in both cases. 

A parenthetical remark is in order. The statement that p = corresponds to outgoing 
boundary conditions should not signify that the resistivity coming from the "bulk" of the 
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sample is directly related to the scattering at the boundaries; the latter are there merely to 
simulate the boundary conditions that enable the "bulk" scattering to produce resistance. In 
other words, the surface should have a roughness profile sufficient to guarantee the boundary 
condition p = 0; however, further increasing the amplitude of the roughness, while increasing 
the resistivity due to the surface scattering, does not increase the resistivity coming from the 
bulk scattering. Therefore, for the case at hand, the actual scattering does not come from 
randomly situated impurities but from the potential of the positively charged background 
ions that form a finite, but otherwise perfect, lattice. A perfect periodic lattice has no 
resistivity, and it is easy to overlook the fact that electrons are scattered by it, i.e., they 
undergo Bragg scattering. As we will show in Sec. [V|, a finite but otherwise perfect lattice 
scatters electrons for all momenta, with two main contributions: constructive interference 
or Bragg scattering as well as background diffusive scattering. It is the latter that leads to 
electrical resistivity and vanishes in the infinite-thickness limit; unlike the case of impurity 
scattering, its resistivity is inversely proportional to the thickness of the film. 

In the remainder of this paper, we will show the details of the calculation of the electrical 
resistivity of a thin film using the Kubo formula within the free-electron model. Specifically, 
we will show that outgoing boundary conditions do indeed imply the existence of a finite 
electrical resistivity. Remarkably, the free-electron result agrees in form and reasonably well 
in magnitude with the one calculated ab initio by Blaas et al. |],|2| for a perfect slab of 
copper embedded in copper. 

III. KUBO FORMALISM FOR THE CONDUCTIVITY OF A THIN METALLIC 

SLAB 

As discussed in Sec. [TI|, a perfectly diffuse boundary amounts to outgoing boundary 
conditions, which are implemented with the corresponding infinite-medium retarded Green's 
functions. As the scattering ultimately leading to resistivity is due to a potential V(r) with 
"finite periodicity," we resolve the Hamiltonian in the form 
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H = H + V(t) , (1) 

where the unperturbed Hamiltonian Hq corresponds to free electrons characterized by 
the eigenf unctions cp^ = fi -1 / 2 e* k r ; in all quantum-mechanical computations, we will use 
particle-in-a-box normalization with finite volume Q. However, given a thin film of cross- 
sectional area A, thickness L, and volume Q = AL, the in-plane dimensions (defining the 
area A) will be effectively regarded as infinite, whereas the finiteness of the perpendicular 
or "longitudinal" dimension L will become the source of finite-size effects; this longitudinal 
direction will be chosen to correspond to the z axis. Accordingly, whenever appropriate, a 
generic vector V will be resolved into its in-plane component Vy and its longitudinal com- 
ponent Vi = zV z ; with this notation, the energy of the free-particle or unperturbed state 
of momentum fik is 

£k = ek||fc z = £k|| + efc z , (2) 

where e k|| = h 2 k\\ 2 /2m and e^ z = h 2 k 2 /2m, with m being the effective electron mass. Our 
calculation of electrical conductivity will be performed by applying perturbation theory 
within the framework of the Kubo formula |7|,[8|]. The required potential matrix elements are 
given by 

= i / d\ e^C"-*)* K(r) = Iv(k- k') , (3) 

where V'(k) is the Fourier transform of the potential. Then, to second order in perturbation 
theory, the diagonal momentum-space elements of the t matrix are 

*k(e) = V kk + £ VWGk,(e)V^ k , (4) 

k' 

where G°(e) is the free-particle propagator. In order to evaluate the Kubo formula for the 
electrical conductivity, the negative imaginary part of the t matrix A = — Im(t) is needed. 
Then, from the condition V^, = Vk'kj it follows that 

A(k h k z ;e) = A(k; e) = -Im [t k (e)] 

^vr^^k-kO^e-ekj-e^), (5) 
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in which the effect of the potential on the conductivity is now summarized by the function 

' V(q) | 2 



^(q) 



Q 2 



(6) 



with q = k — k', and where the density of states S(e — — e^) has been directly extracted 
from -Im[Gl, k ,(e)]/n. 

The Kubo formula |7| gives the zero-temperature in-plane (CIP) dc conductivity in terms 
of the in-plane current-current correlation function, by means of the formal double-limit 
expression M 



a = lim lim 







dre™ <r T j B (r)-j].(0) 



(7) 



which in the independent-electron approximation reduces straightforwardly to the familiar 
form 

h 

: 5{e F - e k - e fe J , (8) 



e 2 1 

— y - 

k|i,fe z ,s 



5e 



d(hh\ 



2A(k h k z ;e F ) 

where s stands for the electron spin index. Equation (§) can be conveniently rewritten by 
evaluating the group velocity as 



dei 



d{Kkn 



2ei 



m 



and applying the in-plane continuum limit 



I v 

if 

kii ,s 



rf 2 k|| 

(2tt)' 



(9) 



(10) 



with the only restriction that the integrand be a spin-independent function of the momentum 
variables |k||| and k z (but not of the corresponding angular in-plane variable). In Eq. (|10|). 
g\\ (ekj is the in-plane two-dimensional (2D) density of states per spin degree of freedom and 
per unit area, which is given by 



011 00 



1 2m 



(11) 



provided that the in-plane dimensions be effectively infinite. Notice that this density of states 
<7ii(e) is actually a constant, g», which from now on will be factored out of the corresponding 
integrals. Then, Eqs. d8|)-(PD imply that the conductivity is given by the expression 
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e 2 1 _ 



h L k z A ( k (tF - e kz ),k z ;e F ) 



(12) 



where k(e) is the positive root of the equation e k = e, i.e., fc(e) = s/2me/%\ in particular, 



M e F - £k z 



1.2 _ L.2 
Hp K z . 



(13) 



In Eq. ([T2]) , it is assumed that A(k||, k z ; e F ) is a function of |k||| and fc 2 alone; thus, to 
simplify the notation, from now on this quantity will be represented as A(|k|||, k z ; e_p). 

Equation (p!2|) will be the starting point for our conductivity calculation in Sec. [V], where 
we will straightforwardly apply its counterpart in the longitudinal continuum limit, 



I roo rile 



1 f°° r 

-y o de kz g ± {e kz ) [/(fc(e*j) +/(-A;(c fc j) 



(14) 



where 



J_/2m\ 1/2 _ 1/2 

2tt U 2 J 



(15) 



is the longitudinal density of states per spin degree of freedom and per unit length; notice 
that, explicitly, k{ek z ) = \k z \. Under this approximation, Eq. (|i~2| ) turns into 

•2 



° = 2h Jo dtkz ^ F ~ tk ^ 9± ^ kz 



1 



+ 



1 



A( k(e F - e fc J, k(e kz );e F ) A( k(e F - e kz ), -k{e kz ); e F ) _ 



(16) 



IV. FINITE PERIODICITY 

A thin film can be regarded as built out of primitive cells assembled into an effectively 
infinite arrangement in two directions (for which we will use symbols 1 and 2) and a finite 
layering in a third direction (for which we will use the symbol 3). As we will consider 
a generic "finite Bravais lattice," the three directions need not be perpendicular to each 
other. In other words, if the number of primitive cells stacked in direction j is Nj, then 
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N±,N2 3> N$, and, in practice, we will regard N\ and N2 as effectively infinite but N% 
as a finite number; notice that iV = NiN 2 N 3 is the total number of primitive cells in the 
film. The finite periodicity in the third direction can be described by counting the number 
M = N% + 1 « N% of stacked infinite lattice planes or monolayers; notice that, for the sake of 
simplicity, we will assume M 3> 1. For instance, if the film has thickness L and consists of 
exactly M monolayers from one boundary to the other, then L = N%d ~ Md, with d being 
the distance between consecutive monolayers. 

The basic periodicity of the crystal structure within its boundaries can be described 
by means of a finite generalization of the concept of Bravais lattice. As usual, given the 
primitive translation vectors a,, with j = 1,2,3, it follows that the set of all translation 
vectors R is of the form R = n\B.\ + 722^2 + n 3 a 3> where the integers rij are limited to the 
values < rij < Nj — 1 for a finite lattice. Then, for any local function of the position, such 
as the lattice potential, the property 



remains valid within the boundaries of the film. In particular, for Fourier transforms, finite 
periodicity guarantees the identity 



where an arbitrary primitive cell C can be related to any other primitive cell via a translation 
Tji by a Bravais lattice vector R. This property leads to the characteristic Bragg peaks with 
respect to electronic conduction, as shown below. In effect, the Fourier transform of the 
potential becomes 



= ^ ^ e - iq - (niai+n2a2+nsa3) / d 3 r e _iq - r y(r) , (19) 

ni=0 712=0 n 3 =0 C 

where V is the entire sample and C is a reference primitive cell. Then, summing the geometric 
progressions involved in Eq. (|19|) and replacing in Eq. (||), one finds 



V(r + R) = V(r) 



(17) 




(18) 




iv 

Ni-l N2-I JV3-I 
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^(q) 



0=1 



^ (0) (q) 



(20) 



where 



^ (0) (q) 



(fir 

— e^Vir) 

c v 



(21) 



is the corresponding quantity for just one primitive cell, v = &x ■ (a 2 x a 3 ) is the volume of 
a primitive cell, and 



sm(N^/2) 
Nj [ sinfo/2) 



(22) 



stands for the interference factor associated with Nj identical primitive cells aligned in the 
direction of the primitive translation vector a 3 , with 



6* = q ■ a, . 



(23) 



The interference factors are the origin of the Bragg scattering peaks, which are represented 
by 5 functions and amount to the selection of the conditions £j = q • B.j = 27mj (with rij 
integer numbers), as can be seen from the vanishing of the denominator in Eq. (P2"|). These 5 
functions can be explicitly displayed by means of the expansion in a series of simple fractions 

1 ^ 1 



E 



sin 2 (e i /2) (6/2 - %tt) 2 

and from the familiar asymptotic result 



(24) 



(0/2) 



(25) 



for Nj — > oo; then, the interference factors become 



sm{N^j/2) 



j oo 

iV 2 ^ 

J n,=— oo 



27T °° 

- £ 5(^-2™,) 

Jv i n,=-oo 



(26) 



where the second expression is to be understood as the asymptotic form for Nj "sufficiently 
large." The continuum limit implicit in Eq. ( |2"6| ) can be reversed by replacing the Dirac 
delta function by a Kronecker delta and recalling Eq. (|23f) , whence 
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^Njitj) ~ E K^^nj ■ (27) 
rij=— oo 

In particular, the combination of the three structure factors in Eq. ( p0|) amounts to 



(28) 

j'=l m,n2,n3=— oo 

Equation (|28"|) can be reinterpreted by expanding q in terms of primitive reciprocal vectors, 
i.e., q = v\b\ + z/ 2 b2 + ^^3, with • b^ = 27c5jh, whence Vj = q ■ aj/27r. Then, Eq. ( |28|) 
states that Uj = rij is an integer, so that q is indeed a reciprocal-lattice vector G nin2n3 , i.e., 

3 00 

n^&)~ e <w m . (29) 

j=l m,n2,n3=— oo 

However, in a finite lattice, it is legitimate to apply the limit of Eq. ( p9|) only with respect 
to the directions defined by the primitive vectors ai and a 2 . Then, the component of q on 
the plane spanned by the primitive reciprocal vectors bi and b 2 is 

(q- a 3), ^ , Qn ^ 

Qin = q 7: b 3 = G nin2 o , (30) 

which is a 2D reciprocal vector, 

gnma = G nina0 = n x bi + n 2 b 2 , (31) 

so that 

2 00 

n^(fc)~ E ^n^n,- (32) 
j'=l n\,Ti2=— oo 

Next we will assume that, for q = k — k', with k and k' on the Fermi surface, the only 
allowed 2D reciprocal-lattice vector g nin2 is the zero vector. To see that this assumption 
is reasonable, let us consider, for example, a thin film of copper, for which one can choose 
the set of primitive vectors ai = (a/2) (x — y), a 2 = (a/2) (x + y), a 3 = (a/2) (x + z), and 
primitive reciprocal vectors bi = (27r/a)(x — y — z), b 2 = (27r/a)(x + y — z), b 3 = (47r/a)z, 
with a m 3.61 A and kp ~ 1.36 Ar x . Then, for q = k — k', with k and k' on the Fermi 
surface, it follows that 
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|q| < 2k F < |bi|, |b 2 



(33) 



implying that the only acceptable choice is n\ = n 2 = 0, namely, 



0. Under these 



conditions, q in = and also qii = 0, because q • &i = and q • a 2 = simultaneously; thus, 
q = q_i_. Then, the argument of the interference factor ^^3(^3) is 



£3 = q± • a 3 = q z d , 

with d being the distance between consecutive monolayers. As a consequence, Eq. 
becomes 



(34) 



(35) 



where 



J i±) (q z )=J r NMd)J {0) (q z 



(36) 



with 



J {0) (q z z) 



d 3 7 



C V 



■V(r)e 



'iq :Z 



(37) 



[cf. Eq. (pT|) 1. Equations (p5|)-(p^) give just a Bragg scattering contribution and zero 
resistivity as N3 — ► 00, but fall short of that singular behavior for N3 finite. Based on the 
preceding analysis, the Bragg scattering contributions, ^ ra gg(q) and J^l gg (q z )-, are defined 
to be the asymptotic forms of the functions J(c\) and J^'{q z ) as N 3 — > 00; thus, they are 
related by Eq. fl3"5]), with 



j£L(q z ) = ^- £ 5( qz d-2nn 3 )J^( 2 ^ 

3 n3=— 00 



V d 



(38) 



in particular, 

i7Bragg(q|| = 0, q z 



T {±) (a 

u Bragg V± 



\q z \<2n/d J\ 3 



(39) 



\q z \<2ir/d 

an expression that will be important in the derivation of the electrical resistivity in the next 
section. 
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V. FINITE ELECTRICAL RESISTIVITY 



Equation fljjp will give a nonzero electrical resistivity only when the t matrix develops a 
nonzero imaginary part, i.e., Aj t ..jfe jt (e) 7^ 0. This imaginary part may arise in the process of 
replacing discrete sums by energy integrals [according to the rule defined by Eq. flljj) ], if the 
self-energy acquires an analytic structure characterized by a branch cut that is effectively 
generated by the merging of the discrete poles of the discrete self-energy. However, the 
scattering by a periodic lattice generates constructive interference in discrete directions 
(represented by S functions); this Bragg scattering fails to produce an imaginary self-energy. 

Due to the directional nature of Bragg scattering, the term J^l gg {k z — k' z ) given in 
Eq. fl38|) does not contribute to the sum of Eq. @ , namely, 

- E ^B r a gg (k-k')Im[G k , ||fc ,(e)] = 0. 

k' h 1 

Then, we are led to a resolution of the function J7"(q) of Eq. @ into two parts, 

= ^Bragg(q) + Jdiff(q) 

= j£U<lz) + ^(Qz)} *q„0, (40) 

corresponding to Bragg scattering, given by Eq. (j38l) , and the remainder, which we identify 
as background diffusive scattering. This procedure, which is based on the analysis of Sec. [IV], 
amounts to isolating the diffusive part of the scattering J7difr(q) = <^"(q) ~~ i7Brag g (q), which 
leads to a finite resistivity via the term 

A(k||, k z ] e) = 7r E ^difr(k - k')S(e - ev - e k >J 

= tE^(*, - KW - c k| - e K ) . (41) 

Equation ([O) can be further simplified by either applying the longitudinal continuum ap- 
proximation, Eq. ( |Hp with respect to z', or explicitly rewriting the 5 function as 

5(e - e k ) = 2irg^\e - e k|| ) \ [5( k z - k(e - e k|| )) + 5( k z + k(e - e k|| )) ] 

**.,k(e-ek„) + ^-fc(e-e k ..) • (42) 
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Then, 



A(k h k z ;e) 



7lL 



-g±{e - e k|l 



J$ ( k z - He - e k|| ) ) + Jtf ( + *(e - e k|| ) ) I . (43) 



Finally, for the calculation of the conductivity, Eqs. (|T|) and (|16|) dictate that Eq. ( f43|) be 
evaluated for electrons on the Fermi surface, for which the conditions \k z \ = k(ep — ekj and 
fcii = k(e F - e kz ) apply; then, 



7lL 



A( A;(e F - e fc J, k z ; e F ) = — <u(e fc J [J^(0) + J$(2k, 
Substitution of Eq. (|§ into Eq. ([16]) leads to 



(44) 



(45) 



where the properties J~^(±2k 



J7difr (2k z ) have been applied. Finally, the term 



J&g{2k z ) can be approximated with its value jj^s (0) = <^diff(0), because it is partially sup- 
pressed by the numerator as k z approaches kp (in fact, the exponential in the corresponding 
Fourier integral does not complete one entire cycle even as z approaches d and k z approaches 
kp and is approximated as taking the value one); then, 



(46) 



h 2nLJ diS (0) ■ 

In order to evaluate the final conductivity expression of Eq. ( f4*6|) , the value of 

J" difr (0) = lim \j (q,| = 0, q z ) - Jsragg (q|| = 0, q z )] (47) 

is required. This can be obtained from the interference factor, Eq. (|22|), which can be 
approximated for N 3 large and £3 small via the expansion 



(£3 



1 sin 2 (iy 3 e3/2) 

m te/2) 2 
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iV 3 



3iV 2 



sin 2 ^ Uofe 2 ), 



(48) 



where the denominator has been expanded in power series of £3/2 and Eq. ( P5] ) has been 
applied to provide the asymptotic form of J 7 ^ (£3)- The derivation above focuses directly 
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on the £3 — > limit and does not emphasize the resolution of the spectrum into infinitely 
many Bragg peaks, each of which has a long tail that contributes to the background diffusive 
scattering. Instead, a more illuminating approach would be to display all the Bragg peaks 
explicitly from the start by considering the limit £3 — > of the simple fraction expansion of 
Eq. ([Zip, namely, 

1 1 



lim 



2 00 1 2 1 

?Sr? M = i' (49) 



,sm 2 (? 3 /2) 2 (6/2)2 

a result that is in agreement with Eqs. ([470 an d (EH)- In the formulas above, in the limit 
M — > 00, the factor sin 2 (M£ 3 /2) oscillates fast about its average value 1/2, which would 
therefore be effectively achieved in all physical measurements; thus, one concludes that 
Eq. ( f48| ) is simplified to 

9tt 1 

^Vafe) = w 6(&) + —j + 0(Q . (50) 

3 3 

Notice the characteristic appearance of the Bernoulli number B 2 = 1/6, associated with 
either the expansion of Eq. ( (48|) or with the value ((2) = B 2 tt 2 of the Riemann ( function 
in Eq. (g^). Then, for £ 3 = q z d small, Eqs. ©, and (0) imply that 

J( qil =0 iqz ) = f N3 (q z d) J®(0) 

= [.Wq,, = 0, qz )} iqzi<2n/d + -^.7< >(0) + Otf) ■ (51) 

therefore, from Eqs. ( p7[) and (|47|), and from iY 3 « M, the diffusive part of J{0) for M 
monolayers is 

^(0) = ^jJ<»»(0) = ^W 2 , (52) 

with 

(v> = / — V(v) (53) 
Jc v 

being the average potential. 

Finally, replacing Eqs. ( fyty and ( |53"D into Eq. PBj), and recalling that L m Md, we get a 
remarkably simple expression for the resistivity, 
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Ti>r 2 d 

p "s7^' (54) 

where 

r = ±J*L vlt)= (m\, (55) 

e F Jc v \ e F I 

which is the average potential relative to the Fermi energy, is a dimensionless parameter 
characterizing the relative strength of the periodic potential. 

Equation ([54]) can be evaluated numerically by introducing a natural atomic resistivity 

p = 4(1 A) ^ 258 /xficm, (56) 
e z 

whence 

270T 2 d[A] 

P~ /ificm. (57) 

Equation (|57D summarizes one of the main results of this paper. It displays a characteristic 
inverse proportionality with respect to the number M of monolayers. In addition, notice 
that: 

(i) The period d is not a free parameter, as it is uniquely determined from the crystal 
structure. 

(ii) As the atomic length scale d[A] is always of the order of unity, the largest variations 
in the resistivity scale of Eq. fl5"7|) will come from the dimensionless parameter T. 

(ii) T is indeed the only free parameter within the framework of approximations used in 
this model. 

(hi) The value of T can be independently estimated from calculations of cohesive energy. 

For example, for copper, reasonable estimates are provided by the following values of 
the relevant parameters §|: d ~ 1.8 A and r ~ 2/7; then, the predicted CIP resistivity is 
approximately 40/M pficm, a value reasonably close to the 116/M /iHcm found by the ab 
initio method [0. 
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VI. CONCLUSIONS 



It is generally recognized that a perfect infinite periodic structure leads to no electrical 
resistance because of the condition of constructive interference or Bragg scattering. In this 
paper, we have shown that, at sufficiently low temperatures, for sufficiently clean samples, 
and for outgoing boundary conditions, finiteness of a metallic film will prevent the potential 
from being perfectly periodic, will produce an effective background diffusive scattering, and 
will cause a size-dependent resistivity inversely proportional to the number of monolayers. 
Our free-electron estimate is in close agreement with similar results from ah initio calcula- 
tions. It would be interesting to see if this resistivity can be measured experimentally by 
properly simulating outgoing (perfectly diffusive) boundary conditions. 

As a coda, it is worth mentioning that the Boltzmann equation approach with the 
conventional relaxation time approximation fails for transport through a finite, yet oth- 
erwise perfect, film. The usual ansatz for the nonequilibrium distribution function 
/(k) = fo(k) + k • E C(k) is not applicable, because the deviation from equilibrium depends 
on the orientation of the electric field relative to the crystal axes of the film as well as on 
the angle between k and E. 
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